Loop Quantum Cosmology of Diagonal Bianchi Type I model: simplifications and 

scaling problems 
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A simplified theory of the diagonal Bianchi type I model coupled with a massless scalar field 
in loop quantum cosmology is constructed according to the p scheme. Kinematical and physical 
sectors of the theory are under good analytical control as well as the scalar constraint operator. 
Although it is possible to compute numerically the nonsingular evolution of the three gravitational 
degrees of freedom, the naive implementation of the p scheme to the diagonal Bianchi type I model 
is problematic. The lack of the full invariance of the theory with respect to the fiducial cell and 
fiducial metric scaling causes serious problems in the semiclassical limit of the theory. Because of 
this behavior it is very difficult to extract reasonable physics from the model. The weaknesses of the 
implementation of the p scheme to the Bianchi I model do not imply limitations of the p scheme in 
the isotropic case. 

PACS numbers: 04.60.Kz, 04.60.Pp, 98.80.Qc 



I. INTRODUCTION 

During the last years loop quantum cosmology [l| 
(LQC) has become an attractive area of research for cos- 
mology and quantum gravity community. The loop quan- 
tum gravity [2|, 13, U] (LQG) which inspired quantization 1 
of the symmetry reduced cosmological models, allows one 
to address fundamental questions about the fate of classi- 
cal singularity and quantum gravitational corrections in 
the early universe. While the topic of the isotropic and 
homogeneous sector of the LQC originated in [6| is well 

understood [IE EE EE EI IMS! there stm is work t0 

be done in the homogeneous, but nonisotropic sector. Al- 
though loop quantum dynamics is not fully understood in 
this sector, already the first calculations in the quantum 
homogeneous models @, H| suggest a completely different 
structure of the space-time near classical singularities. 
After the discovery of the so-called "improved scheme" 
in [l6j , analytical issues of kinematics and dynamics were 
studied in detail in the case of the simplest homogeneous 
diagonal Bianchi type I (BI) model [1J|. Although there 
has been recent progress in the LQC diagonal BI model 
at the level of effective equations of motion [I EE EH , 
the predictions coming directly from the quantum theory 
are still missing. Potential numerical simulations of the 
quantum dynamics in the theory constructed in pjj] are 
more difficult than in the isotropic case [IE] because of 
the complexity of the quantum scalar constraint. While 
the numerical study of the [lj] model is still to be done, 
we can in the meantime focus our attention on recent 
papers [H, EH • In 03 it was shown that a small simpli- 
fication in the quantum theory leads to an exactly sol- 
uble LQC model. Similar simplifications enable us to 
prove self-adjointness of the quantum Hamiltonian in a 
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1 To learn more about deep relation between LQG and LQC see 
introductory review ||. 



purely analytical manner [12j . This motivates us to inves- 
tigate similar simplifications in the LQC diagonal Bianchi 
I model. 

Similar results have also been obtained independently 
by the Madrid group [H, 

Furthermore, so far in the literature there are two kinds 
of " loop regularization" of the gravitational part of the 
scalar constraint for the BI models described in [ll[, the 
so-called fx and fx! . The fx' scheme has better scaling 
properties (see [111]); however, a quantum theory for this 
scheme is much more difficult to construct. In this paper 
we study the fx scheme quantization described in [14| and 
simplify it to some extent. Unfortunately, such a theory 
has a limited domain of applicability. Full scaling invari- 
ance with respect to the different choices of the fiducial 
cell is broken, which causes fiducial cell and fiducial met- 
ric dependence on the semiclassical limit for the fiducial 
cells (and fiducial metric) different from cubical (see Sec. 
IIII Bl for a more detailed discussion). These problems are 
especially important for the case of noncompact topology 
of the three-dimensional spatial slice S of the foliation, 
where different shapes of the fiducial cells (and fiducial 
metrics different than isotropic) are allowed. The quan- 
tum theory should be invariant with respect to changes 
to different fiducial cells (and fiducial metrics) just like 
the classical theory does respect this invariance. The 
model here does not respect the fiducial invariance, re- 
sulting in the semiclassical limit of the theory is being 
unfortunately ill defined. 

This paper is organized as follows. In Sec. II the 
classical theory of BI cosmology is described in terms 
of Ashtekar variables. In Sec. Ill we consider the quan- 
tum theory of BI in the volume and the connection " (77)" 
representations. Section IV describes a unitary transfor- 
mation W which allows one to use continuous 3D Fourier 
transform in order to extract semiclassical states from the 
theory. Section V contains a discussion. In the appendix 
the reader can find numerical strategy used in the simu- 
lations. 
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II. CLASSICAL THEORY 

The scalar constraint for the (minimally coupled) gen- 
eral relativity with a massless scalar field in terms of 
Ashtekar variables is given by 



C' = 



1 



16ttG 



- 2(1 + 7 >- 1 i^A^ ] J + C , (2.1) 

where detE = e 2 and X* stands for extrinsic curvature 
of the spatial slice S. Symplectic structure is defined by 
the following Poisson bracket: 

{A{x)\,E{y)]} = SnG^J^y) . 

Spatial part of the metric tensor for the diagonal Bianchi 
I models is given by 

q ab = a\{t) %l %l+al(t) ^ % 2 b + a 2 3 (t) "bjj %l , (2.2) 

where the left invariant one-forms satisfy <9[ a — (i.e. 

structure constants Cj k =0). If we want to pass to the 
Hamiltonian formulation for our BI model we must be 
careful, because the scalar constraint (12. ip for the metric 
(12. 2p is infinite. In order to make the integration over 
spatial slice finite we can replace £ — > V as 



(2.3) 



where V is a finite cell or, in the case of Bianchi I, we can 
choose topology of £ to be a three-dimensional torus. 
The metric tensor given by the (|2.2p reduced symplec- 
tic structure in terms of Ashtekar variables is defined as 
{c l ,pj} = S^&nG^V^ 1 , where c l ~ 70,^ and pi = ±a,jCLk 2 
for Eijk = 1. Vq = L\L 2 L 3 denotes the fiducial volume of 
the finite cell (or a volume of the 3-torus) . After symme- 
try reduction given by the metric (|2.2j) constraint (|2.1j) 
is reduced to the following one 

C> = - N (t' h I 9 ?, 1 . (cVA + cVic 3 ^ 



87rG 7 2 VbiPaPsI 



+ cWp 3 ) + W(*')o 



1 n 2 



2 



= 



(2.4) 



The variables c z and p 3 are now rescaled as follows 

c 1 = Lie 1 pi = p\L 2 L 3 = ±L 2 a 2 L 3 a 3 . (2.5) 
The symplectic structure is now given by 
{c\ Pj } = 5}8nG 7 . 



Moreover for the scalar field in (|2.4[) we have 
{(j>, n^} = 1. If one chooses the lapse function to be 
-^(£')|piP2.P3| 1//2 = 1 the scalar constraint becomes 



C 



'g^^(cVlC 2 p2+cVlC 3 p 3 +C 2 p2C 3 p 3 ) + ^n 2 = . 

(2.6) 
(2.7) 



Hamiltonian equations of motion 3 

? = {c\C} Pj = { Pj ,C} 

ri = {n ,c} <j> = {<t>,c} 

give us the following conditions: 

c^'pi — constj n«j4 = const^ , (2.8) 

where (z) means no sum over i. Let us set const^ := 
SirG^KICi and const^ :— h^^nGJC^. The solutions can 
be written as functions pi — p(<f)i given by 

-1 - Ki 



Pi ((l>)= PO iexp(VSTTG *-(</>- (foj) , (2.9) 

where Ki = Kni and K^ — Kn^ are rescaled such that 
K X + k 2 + k 3 = 1 (2.10) 

2 1 — i 2 1 2 1 

K 1 + K 2 -f- K 3 -f- Ki — 1. 



There are two types of solutions: " Kasner-like" , when 
two Ki are positive and the other is negative, and 
"Kasner-unlike," when all three K>i £1X6 positive. Let us 
now define, for the purpose of the quantum theory, a new 
symplectic structure as 



{ V \V j } = 12nG^ j , 
where new variables are defined as 



(2.11) 



1 



\/\P(i)\ 



V t =sgn( Pi )H 3 / 2 . 



The classical solutions (12.91) in terms of new variables Vi 
are in the form 

Vitf) = V 0i exp ( ± V8^g| (1 ~ ^ (0 - O )) , (2.12) 

where we put = ±|/e^|, because can be positive or 
negative. The total physical volume of the fiducial cell 
(or a 3-Torus) is defined &s V = IV^V^ 1 / 3 . 



III. QUANTUM THEORY 

A. Kinematics - Volume representation 

A quantum theory in the improved (or not improved, 
so-called °/i) scheme was constructed in detail in [T3 |. 



Plus and minus mean two possible orientations of the spatial 



triads e' 



-a(t) 



3 A time derivative in l|2,7| l is taken with respect to the time t 
associated with laps N(t) = |piP2P3| 1//2 
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In this subsection we briefly recall the kinematics of 
the diagonal Bianchi I model in the so-called ft scheme. 
The kinematical Hilbert space is given by H.Kin = 
L 2 (MlBohr, o^Bohr)® 3 with the orthonormal basis elements 
labeled by three real numbers 

Wi, v%-, va) ■= Wi) ® K) ® \tr 3 ) , 

where RBohr stands for the Bohr compactification of a 
real line. The kinematical scalar product is defined as 



^1,^2, ^3^1, "2, V3) 



, V\ ^V2 , VI ^V3 , V3 ■ 



Any state l^) € H-Kin can be decomposed in the orthonor- 
mal basis as 

Vx V 2 V3 

with the norm N given by 

V\ V2 V3 

There exists a volume operator defined as 

^1,^2,^3) = 67rv / A7Zpi|^ii/ 2 ^3| 1/3 ki,^2,^3), (3.1) 

which is constructed from the components of the volume 
operator as V = | Vi V2 V3 1 x / 3 . The constant A will be 
determined later. Each of the Vi operators is given by 



Vi\Ui,V2,U 3 ) = 67TV AjlpiVi I Vl,P 2 ,U 3 ) (3.2) 

for i = 1,2,3. The operators corresponding to the con- 
nection components Ci do not exist in 7iKin, but there 
are important unitary shift operators (for i — 1,2, 3) 



U$ = exp (ib 



M(i)< 



(3.3) 



where b S R (and there is no sum over i). Their action 
on a basis element is defined as follows 



u£)\vx,v 2 ,vs) =exp {ib^—)\v u v 2 ,v 3 ) = \vi, v 2 + b, v 3 ). 



M2C 



B. Dynamics - Volume representation 

In the paper [l4[ one can find a detailed construction of 
the operator corresponding to the (12. 4p in the full Loop 
Quantum Gravity scheme. The full LQG quantization 
of the scalar constraint (|2.ip gives us quantum correc- 
tions to the term e _1 J5"J5|e' J 'fc by the so-called Thiemann 
trick. It also gives quantum corrections to the term F^ b 
by SU(2) holonomy along a suitable loop (see Eqs. (5.2) 
and (5.4) in [l4| for the case of the diagonal Bianchi I 
model). However, the simplified LQC model proposed in 
[l3| can be interpreted as a quantum model, where only 



the LQG effects coming from curvature 2-from F^ b are 
taken into account. In a similar way we construct in this 
section an operator corresponding to (|2.6p where only a 
curvature 2-form F^ b ~ c^'c^'e^ is loop quantized, so 
in order to simplify quantum theory the Thiemann trick 
is simply ignored. Therefore, we can use a regularized 
version of (121|) (described in HEH) 

given by 



l 



8ttG 7 2 



smfiic sm/Z2C 

Pi = P2 



Mi 



/'2 



sin /iiC 1 sm/i 3 c 3 

-Pl = P3 



Ml 

in 2 - 



M3 



sin /i 2 c 2 sin /i 3 c 3 

-P2 — = P3 



P2 



M3 



, 



where 



Mi 



M2 



M3 



(3.4) 



(3.5) 



A stands for the minimal non-zero eigenvalue of the area 
operator [22| . Let us now investigate the last two expres- 
sions in more detail. The (|3.4|) was obtained in two steps. 
Firstly the fiducial cell Vq = L1L2L3 was introduced in 
order to make the scalar constraint (|2.ip finite. Secondly 
the loop regularization with the condition (|3.5p was used 
to obtain the formula (|3.4|) . Because there is no preferred 
fiducial cell V the theory should not depend on the choice 
of the fiducial cell. In other words the two different cells 
should give the same physical evolution. This indeed is 
the case for the isotropic k = model [16|. However, for 
the quantum model determined by (13. 4p the situation is 
different. It can be easily seen from the (|3.4|) . namely the 
expression sin(Mic 1 ) (the scaling properties for the M2C 2 
and M3 c3 are analogous). Let us notice that 



Mic 



a 1 



IPil 1 / 2 



V a 2 fl 3 

Now, if we rescale our fiducial cell as 

Vq = L1L2L3 — > V ' = L1I1L2I2L3I3 = hhhVa 
the (13. 6p is rescaled according to the (|2.5p as 



(3.6) 



Mi c 



Mic 



(3.7) 



The expression Mi cl i s then not invariant under the 
change of the fiducial cell. The sin(Mic 1 ) is promoted 
to the shift operator by Q3.3p . so it seems that the mag- 
nitude of the volume shift in the operator corresponding 
to the (|3.4p is fiducial scale dependent. This property 
can be summarized as follows: The classical dynamics 
does not depend on the choice of the fiducial cell, while 
the quantum dynamics does. Such behavior is very un- 
fortunate, because it can generate several problems as in 
the old isotropic LQC models. One can argue that for 
the 3-Torus topology where the total coordinate volume 



4 



of the spatial slice E is fixed, the p, condition has correct 
implementation [23, 24]. However, in the case of noncom- 
pact £ topology there is no doubt that (|3.7p is a serious 
defect. So what can we do about it? There are at least 
two possibilities. We can develop a better loop regular- 
ization with correct scaling properties or we can assume 
that the only consistent (partial) solution of the rescaling 
problem is to use only cubical fiducial cells. Then for the 
cubical choice we have k = I (i = 1,2,3) and only for 
the rescaling defined by Vq — ► Z 3 Vo (each a.j — > Zaj) the 
(|3.7p and the volume shift defined by jlid (i — 1,2,3) 
are invariant. In the rest of this paper the cubical shape 
assumption of the fiducial cell is made. Moreover, also 
the fiducial metric °q ab — a uo\ %l +°uj 2 °u% °uil must 
be assumed to be isotropic. 

Our task now is to promote p.4[) to a well-defined, 
symmetric operator. Let us start with a matter term 
(l/2)n|. Because we are interested only in quantum 
gravitational effects the scalar field cf> is quantized as 
usual 

where ip belongs to the total kinematical Hilbert space 
defined as Htot — T^Kin <8> L 2 (M., d<f>). The next step is to 
use kinematical tools described in the previous section to 
get an operator corresponding to the function sin jl^c 1 
defined as follows 



The action of the (|3 . 1 0|) on a basis element is as follows 4 
9 

Cg>i, v 2 ,v z ) = -irGh 2 [{u x + l)(i/ 2 + l)h + 2, u 2 + 2, i/ 3 ) 



- [yi - 


1 \ / I 


1) 


V\ 


-2,i/ 2 + 2, 




- \y\ + 


1 \ / 


1) 


V\ 


+ 2,1/2- 2, 


V3) 


+ (yi - 


1 \ / 


1) 


V\ 


-2,1/2- 2, 




+ {V2 + 


l)\V3 + 


1) 


V\ 


V2 + 2, 1/3 4 


- 2; 


- [l>2 ~ 


1 ^ I, , 1 

l)(l / 3 + 




V\ 


V2 - 2, 1/3 4 


- 2) 


~ {V2 + 


- 


1) 


V\ 


v 2 + 2,1/3- 


-2) 


+ (v 2 - 


i){yz - 


1) 


V\ 


u 2 -2,1/3- 


-2) 


+ {yi + 


l)(^3 + 


1) 


V\ 


+ 2,i/ 2 ,i/ 3 4 


-2) 


- (yi - 


1)(^3 + 


1) 


V\ 


- 2,i/ 2 ,i/ 3 4 


-2) 


- {Ul + 


l)fo - 


1) 


V\ 


+ 2,i/ 2 ,i/ 3 - 


-2) 


+ {v\ - 


1)(^3 ~ 


1) 


V\ 


- 2,i/ 2 ,i/ 3 - 


-2) 



(3.11) 

The Hilbert space preserved by the (|3.If [) operator is 
given by 

-ffpiy = Span^ei + 2n) <g> \e 2 + 2m) ® |e 3 + 2k) : 

n,m,keZ\ (3.12) 

and therefore we call it the physical Hilbert space. The 
parameters £j € [0, 1] . The physical scalar product in 
7Y? hy is defined by the kinematical one by restriction to 
the states ip S 7i? hy . 



sm(fi {k) c k )\i/ k ) = ^(exp(i^ {k) c k )- 

expH^ fe) c fe )) |z/ fe ) = + 2) - K - 2)) (3.8) 

for fc = 1,2,3. Let us now fix the factor ordering am- 
biguity of the (|3.4p operator in a following way. Define 
now an operator 



6 fe K): = -( 



1 /sin/i (fe) c fc T> f sin 



(3.9) 



67T7^ 



PI 



2/ 



((z/fe + l)|z/ fc + 2)-(i/ fe -l)|i/fe-2)) 



and notice the following property: [O it 0^]|i/i, 1/2,1/3) = 
for i =/= j. In terms of (|3.9|) operators we get a symmetric 
(with respect to the kinematical scalar product) gravita- 
tional part of the scalar constraint operator given by 



C gI \l/l, V2,V3) 



8ttG7 

Q.lQ. 3 )\l/i,l/2,l/3) ■ 



O0+O0O, 



(3.10) 



C. Connection representation 

In this section we investigate the quantum theory de- 
scribed in the previous section in (dual) rj 1 connection 
representation. Let us begin with the following Fourier 
transform 7^ hy -> - L 2 (S 1 , d/i 5 i )® 3 

1^)= ^(ei+2n,£2 + 2m,£ 3 +2fc)x 

|e x + 2n) ® |e 2 + 2m) ® |e 3 + 2k) -> (3.13) 
ip(r) 1 ,ri 2 ,ri 3 )= ^ ^(ei + 2n,e 2 + 2m,e 3 + 2fe)x 

„;^(£l+2nhV^(e 2 +2mhV^f(e 3 +2fc)r, 3 



where variables 77 1 € [0,2n/ %/A] 5 6 - From now on we will 
be interested in the most important case when e l = for 



4 An operator 113.111 1 can be obtained also for the laps function 
N(t) = 1 by neglecting the inverse triad corrections and by spe- 
cific choice of the factor ordering in the full quantum constraint. 

5 L 2 (S\d/i s i)® 3 stands for the Hilbert space of square integrable 
functions on a (S 1 ® S 1 ® S 1 ). 

6 The (| 3 - 1 3 P > is a 3 dimensional generalization of the Fourier trans- 
form introduced in [T2II (see eq. (40) and (41) therein). 
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i = 1,2,3. The physical scalar product in the W?_q is 
defined in the obvious fashion: 



(VA) 3 
(2tt) 3 



d 3 ijip(rj)(t)(rf). (3.14) 



/0 "'O ^0 

It is not difficult to check that the basis element 
(rf\2n,2m,2k) = e^ 2n ^ e ^(^W 



is normalized to one with respect to (|3.14p . It is easy 
to derive basic operators in the rf representation from 
the symplectic structure. Let us recall that the classical 
Poisson bracket is given by 



{V k , 



2 S- 
pi e 



. e u) 



A > ( 3 - 15 ) 



where /2(fc)C fc = ^/Ar] k . After canonical quantization we 
get elementary operators. Components of the total vol- 
ume operator are in the form 



for k = 1,2,3. The shift operators 
cation operators 



act as multipli- 



e'^^r; 1 ,)? 2 ,?) 3 ) = e 1 ^ - V(?7\ ^W) . (3.17) 
Let us now consider a part of (13.41) . namely, ^ MC 



If we define the Ok operators as 
.sin(v / A?7' c ) d i 



O k 



drf 



cos(A/A?7 fc ), (3.18) 



then the operators corresponding to sin(/i( i jc l )? i in the 

symmetric factor ordering are defined as 1277rip 1 0fe for 
k = 1,2, 3 7 . The gravitational part of the scalar con- 
straint operator (|3 . 1 0[) can now be written as 



Q r = -187TG7T O1O2 + 2 3 + OiO E 



(3.19) 



The operator (|3.19p is now symmetric with respect to the 
scalar product (|3. 14[) . because the commutator [O,-, Oj] — 
for i 7^ j. Moreover, the spectrum of (|3.19p can be 
determined from the following equation 

c gI ^ {X) (v\v 2 ,v 3 ) = -isttGH 2 ^) i> {X) (v\v 2 ,v 3 )- 

(3.20) 

The eigenvalue has a form w(A) = A1A2 + A1A3 + A2A3 if 
the eigenfunction is given by 



7 Notice that 1277rZp 1 C\- = Ofc i n the volume representation. 



where t/j(\ k )(r) ) are defined by the formulas 
fe lexp(iA fe ln|tan^|) 



(3.21) 



sin v Afy 



for fc = 1,2,3. The spectrum of (|3. 19|) is continuous, be- 
cause its eigenf unctions are normalized to Dirac delta 
(VVamIVVX)) = <^ 3 (^' — Followed by the strategy 

adopted in fbH . [l6| we can now write down the quantum 
constraint equation as 



2h~ 2 c g Mv\v 2 ,v 3 



-c^W,^)- (3.22) 



and interpret the scalar field cj> as time in the quantum 
theory. In order to do this, we can now decompose the 
Eq. (|3.22|) into positive end negative frequencies as 

y/Wmv 1 , V 2 , V 3 , <t>) = ± id^rf.if, 0), (3.23) 



where H = —2h 2 C gr . The y\H\ operator generates 
translations in the quantum "time" (f> therefore we call it 
quantum Hamiltonian. The wave function which satisfies 
the two above equations can be written as 



(3.24) 

where ip(X) is a profile of the wave packet, which will be 
determined later. At the end of this section please notice 
an important property of the above solution. The (|3.24|) 
is symmetric with respect to 7r/vA in each variable rf 
because it shares the same symmetry with (I3.21j) . More- 
over, the states which satisfy (|3.23[) are called physical 
states, because of the Dirac quantization program which 
is applied here. 



IV. THE UNITARY TRANSFORMATION W 

A. Isotropic k = model 

In the recent papers [HI, UH it was shown that the 
differential relation 



Adrj/dy = sin(v Arj) 



(4.1) 



leads to a very simple formula for the gravitational part 
of the scalar constraint operator, namely 



d 2 



in [191, or 



dy 



, potential term 



(4.2) 



(4.3) 
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in [13, where (up to VThKJ factor) y = In | tan 8 . 
The range of vAr? is [0, 7r/vA]> so we can forget about 
taking an absolute value | tan 1 = tan because 
tan > . Notice that in flUE]) sin(VA?y) > as well. 



The interval [0,7r/vA] is mapped to a real line R. We 
will see however that in the Diagonal Bianchi I model the 
situation is different. 



B. Diagonal Bianchi I model 

A simple generalization of (|4.1[) for the quantum 
Bianchi I model yields the following differential relation: 



y/Adrj k /dy k = sin(\/A77 fe ), 
which gives us the formulas 



hi 



tan ■ 



Arj k 



(4.4) 



(4.5) 



From (|3 . 13[) however we will find that the range of the 
variables rj k is now [0,2tt/\/A] for k = 1,2,3 ( in contrast 
to the isotropic case where the range is [0,7r/\/A]). Let 
us consider this relation in more detail. Notice that: (1) 
y k (r] k ) is not a bijection; (2) y k {r] k ) is symmetric with re- 
spect to n/VK. Moreover, from (|4.5|) it is clear that the 
interval [0,7r/\/A] is mapped to a real line (denote it by 
R + ) because y k (0) = -oo and y k (n/y/A) = oo. The 
second interval [7r/\/A, 2tt/\/A] is mapped to another 
real line (denote it by R_) because y k {-K/\fK) — oo and 
y fe (27r/v / A) = — oo. Now, because the right-hand side 
of (|4.4p as well as tan m Eq. (|4.5[) are positive for 

■q k £ [0, 7r/VA] and negative for r/ k S [7r/\/A, 2-K/y/A] 
(for k = 1,2,3), it is natural to consider y k (f] k ) in each 
of the two intervals separately. Let us now consider the 
physical scalar product (|3.14[) and the physical states 
(|3.24p in more detail. Because of the very important 
symmetry of (|3.24[) with respect to n/ -s/A in each of the 
variables rf , the physical scalar product (|3 . 14[) can be 
reduced to 



=8 ^f /'"/""/ '"' /:i '7'-('7"'('7) N-«i) 
t^) ./o Jo Jo 



Keeping this property in mind we will now construct a 
quantum theory in the y k variables in the symmetric sec- 
tor for r\ k £ [0, 7r/\/A]. Equations (|4.4 [ I4.5|) define a uni- 
tary map W (which is a three-dimensional generalization 



of the map introduced in [12|]) as 

y ay L y dy z y ay 3 



x V^V^V^V)) (4.7) 



and W as 



1 / dy 1 / dy 2 / dy 3 
(VA) 3/2 V ^ V ^ V 
></(iW),tfV),tfV))- (4.8) 



Let us denote by £(1/2) a half of the circle, namely the 

interval [0, tt/vA]. The map transforms then the 
Hilbert space L 2 (Sf 1/2) , dfi s i)^ 3 to the L 2 (M 3 , d 3 y). The 
scalar product in the new Hilbert space is defined as 



oo poo 



{m = wLLL d * ,d * v 



(4.9) 



Let us now consider the operators. The Ok [defined by 
(|3.18p ] are transformed under W to well-known operators 
as follows: 



W(-i 



sin(v / A?7 fc ) d 



■ cos(VAri k ))W' 1 = -i-Q^j: 



/A di] k 2 

(4.10) 

for rj k £ [0,7r/vA] (where k — 1,2,3) only. The com- 
ponents of the volume operator (|3.16[) are transformed 
under W as 



(4.11) 



W V k W- 1 = 12tt 1 1 2 1 W( - i-^W- 1 = 

I2^l 2 ,v/A( - i cosh(y fe ) A _ i S mh(y k )) . 

Moreover it is not difficult to find a form of any physical 
state p. 241) in new variables 



Our variable \Z~Arj is denoted Xb in [T3l and x/2 in [l2ll . 



iw(A)>0 

(4.12) 

where A • y — Xiy 1 + \iy 2 + A3?/ 3 . 

At the end of this subsection please notice the simple 
fact that if the physical states (|3.24j) share the following 
property ^(ff) =0 at the boundary of the region of the 
integration in (|4.6p then the operators (|3.16|) and (|4.11|) 
are symmetric on L 2 {S^^y dn S i)® 3 and L 2 (R 3 ,d 3 y) re- 
spectively. 



C. Classical and quantum evolution 

In this section we present the differences between clas- 
sical and quantum evolution of the scale factors in the 
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diagonal Bianchi type I model. The classical relation be- 
tween elementary variables and the scale factors [in the 
metric tensor (|2.2|) ] is 



\pi\ = a 2 a 3 
\p 2 \ = aia 3 
\Pa\ = aia 2 



(4.13) 



New variables were introduced, namely Vj = 
sqn(pi)|pi| 3 / 2 , and these have corresponding opera- 
tors in the quantum theory. The classical trajectories 
are defined as 



3(1 
2 l 



K4\ 



(4.14) 



Moreover, the parameters Ki and k^ satisfy relations 
()2.10p . In order to compare the classical and quantum 
model we can use the semiclassical states which are ana- 
logues of those introduced in [l4[ 9 10 , namely 



d 3 X^(X) %)(j7)e = 



*V5 



j(A)>0 

where the amplitude is given by 



/8irG^/Z(4> -<t>*) 



(4.15) 



the non linear behavior of the |w(A)|( 1//2 '. The func- 
tion e lX -y ^(X)e pi ^ V ^° y/ ^ {4 ' ^ ) , however, is under 
good analytical control and we can compute the inte- 
gral (14.18P numerically using fast Fourier transform in 
three-dimensions (see appendix for details). 

From the result of numerical simulations one can no- 
tice an important property of (|4. 18|) and as a consequence 
of the W map, the property of (|4.15|) . Numerically com- 
puted W(9)(y, (f>o) has a finite norm 11 in the sense of the 
scalar product (|4.9j) . Moreover, W(9)(y, </>) decreases as 
cxp(— ayl) (where a > is some constant) in every di- 
rection y k for all values of (f> used in the numerical simu- 
lations. If we use the inverse transformation W~ x 
and write the state (|4.18|) as 



1 



1 



*(rf, 4>) = - _ _ _ 

y'sin(%/A7 ? 1 )) y'sin(VAr7 2 ) yj 'sin(\/zV) 

xW{9)(y(rf),cj>), (4.19) 

we conclude that 9(rf,4>) — at the boundary of the 
integration region in (|4.6[) . In order to make it more 
clear, let us write the right-hand side of (|4.19[) in terms 
of y k variables as 



_ _ (A!-K«i) 2 (A 2 -K«2) 2 (A3-C3) 2 

ip(X) = e < 2 -i e < 2 -2)2 e ^"^' 2 . (4.16) 



Parameters Ki satisfy relations (|2.10p . Consider now 
(|4.15p . which is sharply-peaked at some JCni, for i = 
1,2,3. Since (|4.16p is nonzero only in the small region 
close to A = (/C«i, /Ca«2, ICk 3 ) and moreover, if K,Ki 1 
then 2(AiA 2 + A1A3 + A2A3) w K. 2 k^ > and we can 
approximate the integral (|4.15p by 



9 (if, 



^GJ\uj\(4, -^) 



d 3 X^(X) ^ Cx) (rf)e Tl - 

(4.17) 

because then the variables Xi (i — 1, 2, 3) for which co < 
do not contribute. For the symmetric sector defined by 
(14. 6p the physical state (|4.17p in the yk variables has a 
well known form of Fourier transform in 3D 



d 3 X ^(Xy^e^vT 



&TrGJ\u)\{<t> Q -4>*} 



(4.18) 

with the amplitude (|4"T6|) . The integral (|4T8|) is 
not easy to compute analytically [apart from the case 
when (6n — (/>*) = 0]. The difficulties come from 

9 In [1411 the scalar constraint operator is defined by a different 
factor ordering in the volume representation, so the factor (cpo — 
<fi*^ in 114.151 ) is multiplied by a different constant. See (8.41) in 

10 The state 114.151 1 is not normalized to one. 



RHS = x/cosh y 1 x/cosh y 2 a/cos1i y 3 W(#) (y , 4>) . (4.20) 

When 77^ — > from (|4. 5|) we get y k — > —00, so the right- 
hand side and ^(77, </>) — > 0. Similarly, when ry' — > 7t/a/A 



we get y k —> 00 and again ^(r/, 



0, because 







lim y cosh j/ 1 \/ cosh y 2 a/ cosh y 3 I4 7 (\l/) (y, < 

y fc — >±oo 

(4.21) 

for fc = 1,2,3. This property is very nice, because then 
the operators (|3.16p are symmetric on semiclassical states 
with respect to (|4.6p . Once the semiclassical states are 
known it is possible to compute the expectation values of 
the components of the volume operator (|3.16p on semi- 
classical states, namely, 



(W(V)\W VkW- 1 W(9)U 
(W(9)\W(9)) 



(4.22) 



The scalar product and the components of the volume 
operator are given by (|4.9p and (|4.1ip respectively. The 
above integrals were computed numerically (see appendix 
for details). 

A norm of 114.181 1 can be also computed analytically. 



DISCUSSION ted as functions of the quantum time <f>. In order to 



In the Figs. JTJ, © and (J3} expectation values (|4T22|) 
computed in different semiclassical states (|4.18[) are plot- 
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6 *J nG (0- O ) 



FIG. 1: Isotropic solution for fti = K2 = K3 = 1/3, = \/2/\/3. On the vertical axis Vi denotes the expectation value 
(Vi)(4>) . The dashed lines denote classical trajectories. 



:ooo 

1500 
10O0 
500 





FIG. 2: Kasner-unlike solution for «i = 1/4, ^2 = 1/6, = 7/12, and « 0.7546. V» stands for the expectation values 
(Vi}(<fi). The dashed lines denote classical trajectories. 






FIG. 3: Kasner-like solution for Ki = —1/5, K2 = 3/5, K3 = 3/5, and K<f, = \/6/5. Expectation values (Vi}(0) := Vi are 
compared with the classical trajectories (dashed lines). 



make the pictures clear, each figure contains quantum 
and classical behavior of one degree of freedom V$. Three 
solutions: isotropic, Kasner-like, and Kasner-unlike are 
compared with classical trajectories. In each case far 
away from classical singularity quantum trajectory fol- 
lows closely the classical one up to the region where 
quantum gravitational effects become dominant. Instead 
of following the classical singular solution we have "big 
bounce" to another classical trajectory defined by (|4.14|) . 
While it seems that classical singularity is avoided we 
will see later that unfortunately this result as well as the 
semi-classical limit of the theory does not have invariant 
meaning with respect to fiducial cell and fiducial metric. 
However, having now all types of numerical semiclassi- 
cal solutions of (|4.22|) we can make a phenomenological 
observation. Each curve in Figs. (|1I2I3[) can be well ap- 



proximated by the following hyperbolic function: 

:3 A 



( Vi )(<j)) = V Ot cosh(V8TrG-(—-^)(0-^)) , (5.1) 

Z \K(fi\ 

where values Vbi and <f>° depend on semiclassical states ip 
which are used in computations. Now, we consider nu- 
merical parameters used in calculations and dispersions 
defined by 



4^) 



Typical values of JCki in the Fourier amplitude (|4. 16|) 
used in numerical calculations are as follows: K,Ki fa 
100, 150, 200, 250, 350, 450. Typical behavior of disper- 
sions is plotted in Fig. j4| for isotropic and Kasner- 
unlike solutions. As can be seen, before and after big 
bounce relative dispersions AVi/(Vi) are asymptotically 
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constant. Typical asymptotic values are AVi/(Vi) ~ 0.2 
in our numerical calculations. The integral (|4.22[) and 



V, 




FIG. 4: The behavior of relative dispersion is plotted as a 
function of 4> (This plot stands for isotropic solution, however 
for the Kasner-like case the plot is practically the same for all 
Vi.). 

dispersions were computed using lattices defined by pa- 
rameters Ni = N 2 = N 3 = 256,320,512 [see $£7}}. In 
principle one can use better parameters 2Vj in (|7.7[) and 
higher values of /Ck^ to increase the accuracy and semi- 
classicality of the calculations. However, the require- 
ments for the calculating machine are beyond the param- 
eters of the " SOWA" cluster in the University of Warsaw 
physics department which the author used and whose ca- 
pability allows for maximum values of Ni, N%, and N3 
in (|7.7I7.10|) equations to be 512. Let us now mention 
the difficulties with the Kasner-like solution. Computa- 
tions of small dispersions for this case are unfortunately 
impossible to perform for the author. This is due to the 
fact that in order to have a good approximation in (|4.15D 

/- -/ 

Ju(A)>0 JR 3 

defined by 2(AiA 2 + A1A3 + A 2 A 3 ) « K?k\ > for suffi- 
ciently large Oi parameters, we have to take large num- 
bers K,Ki, which is not possible according to the random 
access memory (RAM) capability of the "SOWA" clus- 
ter. Our maximal lattice defined by Ni = 512 (i = 1, 2, 3) 
is not dense enough to compute integral (|4.22p with very 
large values of /Ck^. However, we do not see a reason why 
dispersions for the Kasner-like solution would be differ- 
ent than in the Kasner-unlike solution, so we expect very 
similar behavior. 

Results enclosed in this paper show only semiclassical 
behavior of the three gravitational degrees of freedom. 
However, the above results are valid only for the cubical 
shape of the fiducial cell and metric. Moreover only the 
cubical scaling of the theory is allowed. Let us now focus 
on the scaling properties of the gravitational part of the 
scalar constraint. If we allow one to scale — > hcii, 
where all U are different for i = 1,2,3 (the fiducial cell 
becomes cuboid) then the spectrum of the operator (|3.10p 
scales as 

a; ~ — 2 (aiAiA 2 + a 2 AiA 3 + a 3 A 2 A 3 ) (5.2) 



where Vq = hhh and ai,o; 2 ,a 3 are some numbers de- 
termined by li, I2, l 3 - For ij ^ 1 we have a, ^ 1 and then 
semiclassical states (|4.15|) do not give correct classical 
limit. The quantum trajectory defined by (|4.22p simply 
does not follow the classical one as one can see in the Fig. 

m 




FIG. 5: The example of the isotropic solutions (Vfe) (for k — 
2, 3) plotted as a function of (j>. The fiducial cell is rescaled 
in one direction as Vo = hhh = ' 3 — > Vo = \ihhh = Ai? 3 , 
where Ai = 2. 

Thus the semiclassical limit with respect to (|4.15|) is 
fiducial cell dependent. Moreover, the shift operators 
defined by sin/XjC 1 correspond to the different steps af- 
ter rescaling Vq — ► Vohhh, m contrast to the k = 
model jTfl, where the step defined by the sin/Zc is in- 
variant under Vq — > 1 3 Vq. One can conclude that for 
each choice of the (noncubical) fiducial cell V (and fidu- 
cial metric) the model has different quantum dynamics, 
while it seems reasonable to expect that the quantum 
model should generate only one physical dynamics. The 
only consistent fiducial cell is a cube with the rescaling 
of the scale factors eij — * Zeij. Then ot\ = cv 2 = a 3 = Vq, 
the spectrum is invariant uj — > to, and the semiclassical 
limit defined by (14. 15)1 is correct. 

Let us conclude that the p, quantization scheme which 
is highly successful in isotropic models, applied for the 
Bianchi I in [lj| has its problems described above. The 
classical theory allows one to consider many different 
fiducial cells, while the quantum model allows only the 
cubical one. Moreover, the fiducial metric 

o„ 0. ,1 0. ,1 1 0. ,2 0. ,2 , a ,3 q ,3 

also must be isotropic, which is a serious limitation. Al- 
though for the 3-Torus topology of the spatial slice S 
with fixed coordinate volume Vq the ft scheme can make 
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sense [23J, l24| , for the noncompact £ topology the model 
is tentative and should be considered as such. 

Now what about future directions for the given model? 
It seems necessary to carry out the analysis involving the 
curvature invariants in order to investigate the scaling 
properties of the model, especially at the bounce. How- 
ever, our results suggest that the £L scheme applied to the 
BI model does not have correct scaling properties even at 
the bounce. We conclude that the new way of loop reg- 
ularization and its quantization is needed. Also similar 
problems may arise in a large class of models based on 
arbitrary discretizations such as lattice refinement mod- 
els, so one should be more careful with respect to the 
regularization of the constraints. 
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VII. APPENDIX 



A. Numerical Integration 



The sum X k := Y^f^ ^. 



Consider a Fourier integral 



dX ip(X)e 



(7.1) 



where the amplitude is sharply peaked at some value of 
Ao and zero everywhere else. In order to deal with the 
above integral numerically we can use the following strat- 
egy. Let us approximate (|7.ip by 



ip(x) 



(7.2) 



If the interval [a, b] covers the region, where the ip(X) is 
nonzero, and if this interval is large enough, then the (|7.2p 
is a good approximation of the integral (|7.1[) . The next 
step is to approximate the integral by defining a lattice on 
the x axis Xk = xo+kA x and on A axis Xj = a+jA, where 
k,j = 0, 1,2... A - 1. Moreover, A x = 2n/(b - a) and 
A = (b — a)/N . The integral can be then approximated 
by the sum as 



N-l 



i){x k ) « A MX 3 )e^ Xk 



(7.3) 



1=0 



N-l 



Ae* ax °e lakA * e l ^ X0 ^{Xj)e l 

3=0 



' , where ipj 



a ijAx 



can be computed using the powerful fast Fourier trans- 
form. As a result we get values of ip(x) on a lattice de- 
fined by x k = x + kA x , for k = 0, 1,2,3. ..N - 1. For 
sufficiently large A (and small A , A x ) it is possible to 
cover the region on the x axis where the function ip{x) 
is nonzero with a good accuracy. A norm of (|7.1[) can be 
computed as 



I = 



dx tp(x)ip(x) 



dx ip(x)ij}(x) 



x 



N-l 



A x 22 ^(xk)ip(xk) ■ 



(7.4) 



k=0 



It is now not difficult to develop a three-dimensional ana- 
logue. Let us consider the following 3D Fourier transform 



oo roc />oo 



i>{y\y\y A )= dX l dX 2 dX 3 

J — oo « — oo J — oo 

x $(Ai, A 2 , X 3 y {Xiyl+X2y2+X3y3) . (7.5) 

If again the Fourier profile f/>(Ai, A2, A3) is sharply dis- 
tributed, the integral can be approximated by 

phi pb2 rbs 

^(y\y 2 ,y 3 )~ dX 1 dX 2 dX 3 

J Cli J 0,2 «/ 0,3 

x ^(Ai, A 2 , A3)e i(A ^ 1+A ^ 2+A ^ 3) , (7.6) 

if the boundary of the integration region defined by <2j 
and bi (for i — 1,2,3) is appropriately chosen. The next 
step is to approximate the above integral by the sum [at 
some point (y^, y% 2 , yfj] as 

Ni—l N 2 -l N3-I 

VG/LWfa.lO^iAaAa £ £ E ( 7J ) 

jl=0 J2=0 j 3 =0 

x ^(xi 1 , x{ 2 , Af y (x ? < + A ^ 2 < +A 3 3 < ) . 

If we set a lattice: 1) y\ x =Vq + hA yi , yf 2 =yl + k 2 A V2 
and y| 3 = y$ + k 3 A y3 ; 2) Af = Oi+Ajji, Af = a 2 + A 2 j 2 
and A^ 3 = a 3 + A 3 j 3 , where Ai = (61 - ai)/Nx, A 2 = 
(62 - a 2 )/N 2 , A 3 = (63 - a 3 )/N 3 , A yi = 2tt/(6i - oi), 



A y2 = 2n/{b 2 - a 2 ) and A y3 = 27r/(6 3 - o 3 ) then 
has the following form: 

Mvk,vL>vk) « A!A 2 a 3 e™^ k ^x klMM , (7.8 

where arg(fci, k 2 , k 3 ) = 2/o°i + y$a 2 + y%a 3 + A Vl kiai -f 
A y2 k 2 a 2 + A V3 k 3 a 3 and 



ATi-l JVa-l JV3-I 

X kM3 = E E E W^Af) 

J"l=0 J 2 =0 j 3 =0 

31 A x + J2 A 2 + y 3 , j 3 A 3 ) 



(7.9) 



x e 



x exp 



.., . k 2 j 2 k 3 j 3 , 
i2-k\— V — h —^-) 



A'; 
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where fej = 0, 1,2..JV< — 1, for i — 1,2,3. Similarly to 
the one-dimensional case the sum (|7.9p can be computed 
using fast Fourier transform in 3D. The scalar product 
(|4.9p is defined as 

g N1-IN2-IN3-1 

m) M 72^w A ^ A ^ A ^ E E E 

V ' fei=0 k 2 =0 k 3 =0 

x^yl^yl^yk^iyk^k^yV ■ (7.10) 
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